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ABSTRACT 


Astronomical observations of SN1987A, such as the light curve, spectral intensities of 
lines, the X-ray emission etc., constrain the lifetime for the decay of a heavy neutrino 
lMeV ~ m„ H < 50 MeV through vg —* i/i-|-e + -|-e~ exceeds 4*10 15 exp(— m Utt /5MeV) 
seconds. Otherwise, resulting ionization energy deposits and stronger x-ray emission 
would have been observed. This coupled with traditional cosmological considerations 
argues that the lifetime of r-neutrinos probably exceeds the age of the universe. This 
in turn would imply the standard cosmological mass bound does apply to v T , namely 
m Vr ~ 100h 2 eV (where h is the Hubble constant in units of 100 km/sec/mpc). The 
only significant loophole for these latter arguments would be if u T primarily decays 
rapidly into particles having very weak interactions. 
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The standard model of weak interactions demands that individual lepton numbers 
be conserved. On the other hand almost every extention of the model envisages 
the possibility of the nonconservation of lepton numbers, thereby allowing neutrinos 
of one type to decay into others. It has been recognized, for more than a decade 
now, that astrophysical and cosmological considerations can provide rather stringent 
bounds on the decay rates of neutrinos through various channels^ 2 ’ 3 ). Among these 
the most stringently constrained mode is the radiative one: i/jy — ► 1/1 + 7. However, 
theoretically, this channel is expected to be suppressed 4 . Alternative modes of decay 
are those with scalars, pseudo scalars, Nambu-Goldstone bosons or e + e~ pairs in the 
final stated. Since the decay products are directly observable the last of these modes 
is of particular interest, especially for the r neutrino whose mass is constrained by 
laboratory measurements to be m„ T < 35 MeV, thus allowing in principle, for above 
the threshold decay into e + e - pairs. The supernova 1987A which exploded in the 
LMC provides a sensitive and direct way of looking for decays of neutrinos without 
the uncertainties caused by assumptions about the cosmological density parameter 
and the supernova rate which go into the earlier discussions. Previous examinations 
(cf. ref. 5,6) have focused on lifetime limits from the pure radiative mode implied 
by the lack of observed 7-coincidences with the */-burst. In this letter we derive 
the constraints placed by the observations in the vicinity of SN1987A on the decay 
VH -» + e + + e~. 

The idea on which these discussions are based is that the electrons and positrons 
that emerge from the decay will lead to observable effects through their interaction 
with the supernova debris and the circumstellar material. These effects are listed 
below: 

1. Deposition of energy in the material through coulomb interactions. 

2. Ionization. 

3. Generation of X-rays through bremsstrahlung. 

4. Generation of X-rays through annihilation of the positrons. 

5. Acceleration of the material by the gradients in the pressure of electrons and 
positrons. 

Limits on neutrinos come from the fact that in gravitational collapse the neutron star 
binding energy, ~ 3 x 10 53 er^s, is radiated in neutrinos. Via neutral currents, all 
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species of neutrinos are radiated. The arguments below will apply to all neutrinos 
with m„ ~ 1.1 MeV so that e + e - pairs can be produced. 

In calculating the bounds we make several simplifying assumptions; more exact 
calculations will not qualitatively alter the limits presented here. These assumptions 
are the following: 

1. The total energy radiated per neutrino species (neutrino plus antineutrino) is 
independent of its flavor 7 and is given by 

Q ~ 10 53 er< 7 s « 6 x 10 68 AfeF for m v = kT v ~ 5 MeV 
* 10 63 (^)e*p-(^) for m I/ »AT 1/ 

[A temperature of 5 MeV is a conservative estimate for /z and r neutrinos. 
Mayle etal. 7 find kT V)t ~ 6 MeV. Only v e and P* have lower temperatures, 
due to their charged current interactions altering the radii of their respective 
neutrino spheres 7 .] 

2. The kinematics of the three body decay are not included in the calculations but 
it is assumed that about a third of the energy of the neutrino is carried away 
by each of the decay products. 

3. We assume that the electrons and positrons that arise in the decay have negli- 
gible spatial diffusion. In the Large Magellanic Cloud, with its high density of 
interstellar gas, the magnetic fields are expected to be somewhat higher than 
the Zfi Gauss fields of our Galaxy, say ~ 5 fi Gauss; the circumstellar fields may 
even be larger. In any case the gyroradii of the electrons and positrons with 
typical energies of ~ 5 MeV in such fields is ~ 10 10 cm. Since the scattering 
meanfree path of particles in a magnetised plasma is about the same as their 
gyroradii the spatial diffusion can be neglected over the length scales of interest 
here. 

With these assumptions we now proceed to calculate the effects. The energy 
radiated in neutrinos is given by 

Q = 10 53 /(m) erg 
/(m) = 1 for m = kT 

/(m) = ( m/kT)exp — ( m/kT ) for m » kT. 


( 1 ) 
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The number of neutrinos emitted in the supernova is given by 

^ - 4 x 10 5 7 e~ m/fcr 


N 


< E v > 


( 2 ) 


where < E„ >~ 3 kT for m ~ kT ~ 6 MeV and < E v >~ m for m kT . The 
number of positrons or electrons generated per unit volume at a distance r is given 
by 


n(r) ~ 


N 


N 


N 


( 3 ) 


4-nrr 2 7/3cr 4irr 2 cr kTAirr 2 cr 
where r is the lifetime for the neutrino decay. Since we are considering neutrinos of 
mass above the decay-threshold of ~ 1.1 MeV , setting 7 « 1 and (3 « 1 in eq. 3 does 
not lead to serious error in our estimates. 

Each electron and positron arising from the decay will have an energy of several 
MeV in the laboratory frame and will deposit energy in the surrounding material at 
the rate of ( « 4 MeV g~ l cm 2 « 4.10“® erg g ~ 1 cm 2 . This will heat up the debris 
and will contribute to the luminosity of the supernova in a manner quite similar 
to the energy deposited there by 58 Co decay, for example. The contribution to the 
luminosity due to this process is given by 

4 • i() 53 e“ m /* T 


L = j i\ = j ^(f, ir = 


3kTr 


-erg s 


-1 


( 4 ) 


Here r 0 « 10 gcm~ 2 is the optical depth of the matter surrounding the supernova 
at present. Now, the present luminosity of the supernova is ~ 10 39 5 er <7 a -1 and the 
light-curve fits a 66 Co decay lifetime of 114 days excellently Study of the light 
curve shows that an additional contribution from any hypothetical process cannot 
exceed ~ 2.5% or = 10 38 er<j a -1 . Substituting this in eq. (4) we get 


r > 4 • 10 15 e~ m/ * r s 


( 5 ) 


Note that even for m„ > kT v ~ 5 MeV this still yields a restrictive limit. 

A sizable fraction of the energy transferred to the debris through Coulomb inter- 
actions of the electrons and positrons goes into ionization of the hydrogenic material. 
At larger distances (10 16 - 5-17 cm) the particle densities are sufficiently low and the 
velocities are sufficiently large, that Lyman-a photons escape relatively freely from 
the material. Under these circumstances a fraction of about 30 % of the ionization 
energy is converted into Lyman-a photons. Thus the luminosity in Ly-a is about a 
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third of the energy transfer rate given in eq. 4 




2 Qir 0 10 ^ e -m/kT 


erg a 


-l 


( 6 ) 


9 kTr t 

Noting that a luminosity of ~ 2.10 38 in Lyman-a would have been detected, even if 
interstellar absorption is taken into account, one gets as before the limit 


r > 10 15 e“ m/ * r a 


( 7 ) 


The limits obtained from the rest of the considerations 3-5 are less stringent; we 
therefore discuss them only briefly. The X-ray luminosity of the supernova is less 
than ~ 10 39 er<; a -1 now( 8 L The expected X-ray flux is estimated simply by assuming 
that the electrons and positrons radiate e -1 of their energy within a radiation length 
A ro d « 100 gcm~ 2 . 

Lx-ray « [^/(m ) ^ 10 39 ers a- 1 (8) 

TArad 

The inequality (8) translates into r > 10 13 f(m) s. The gamma-ray data has been 
obtained with telescopes with very large angular response and yield far less stringent 
limits. The dynamical effects of the electrons and positrons on the debris could be 
used to preclude very short decay times for the massive neutrinos by noting that the 
kinetic energy in the debris is a small fraction of the energy emitted in neutrinos. 
The dynamical effects for relatively large r are subtle and need further study. 

It is interesting to consider the implications of the limits derived here. Among 
these the most stringent is the one derived from the observations of the light curve in 
eq. 5. Let us focus on discussions on the r-neutrinos whose mass is constrained by 
laboratory studies 9 to be less than 35 MeV and thus may be in a regime where these 
limits are relevant. Depending on the actual value of its mass the lifetime for the 
decay into electrons and positrons the r must satisfy the following set of conditions 


r > 4 • 10 15 a for 1.1 < 

m 

w 5 MeV I 

r > 4 * 10 15 a for 

m 

« 10 MeV > 

r > 4 • 10 13 a for 

m 

« 35 MeV 1 


This means that the r-neutrinos generated in the big bang will survive without sub- 
stantial decay up to redshifts (1 + zf) « (T u /r) 3 / 3 i.e. up to redshifts of 100 to 
1000 depending upon their mass. That is the decay will occur after the decoupling 
of the cosmological background radiation. And, if their decay lifetime should be as 


short as indicated by the limits then there would be several important cosmological 
consequences 3 , distortion of the relic microwave background and the effect on the 
cosmological expansion, to name only two of them. The latter effect is discussed in 
detail by Dicus, Kolb and Teplitz who derive the constraint 

m(r/r u ) 1/3 < 100 h 2 eV (10) 

From equation (10) one can see that either the mass of the r-neutrino is below the 
1.1 MeV threshold for t/e ± decay or that its lifetime is longer than the age of the 
universe. Keeping in mind the theoretical and experimental constraints on radiative 
decays this long lived option would mean that the r-neutrino is also very light and 
should satisfy the constraint 

E m Vi < 100 h 3 eV (11) 

unless, of course, it can decay rapidly into particles having only very weak interactions 
with matter. 
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